The importance of electron-electron interactions in the RKKY coupling in graphene 
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We show that the carrier-mediated exchange interaction, the so-called RKKY coupling, between 
two magnetic impurity moments in graphene is significantly modified in the presence of electron- 
electron interactions. Using the mean-field approximation of the Hubbard-t/ model we show that 
the (1-f cos(2ku • R)-oscillations present in the bulk for non-interacting electrons disappear and the 
power-law decay become more long ranged with increasing electron interactions. In zigzag graphene 
nanoribbons the effects are even larger with any finite U rendering the long-distance RKKY coupling 
distance independent. Comparing our mean-field results with first-principles results we also extract a 
surprisingly large value of U indicating that graphene is very close to an antiferromagnetic instability. 

PACS numbers: 75.20.Hr, 75.75.-c, 73.20.-r 
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Several novel features of graphene, such as two- 
dimensionality, linear energy dispersion, a tunable chem- 
ical potential by gate voltage, and a high mobility have 
helped raising the expectation of graphene being a se- 
rious post-silicon era candidate [l|-l3|. In this context, 
functionalization of graphene, especially with magnetic 
atoms or defects which also opens the door to spintron- 
ics is of large interest. One of the most impor- 
tant properties of magnetic impurities is their effective 
interaction propagated by the conduction electrons in 
the host, the so-called Ruderman-Kittel-Kasuya-Yoshida 
(RKKY) coupling 0. This coupling is crucial for mag- 
netic ordering of impurities but also offers access to the 
intrinsic magnetic properties of the host. Several studies 
exist for the RKKY coupling in graphene, where both 
the standard perturbative approach applied to a con- 
tinuum field-theoretic description of graphene and 
exact diagonalization have been shown to give simi- 
lar results. However, consistently, the RKKY coupling 
in graphene has been calculated for non-interacting elec- 
trons. This is in spite of growing evidence for the impor- 
tance of electron-electron interactions in graphene with 
theoretical results pointing to intrinsic graphene being 
close to a Mott insulating state l3"12|. These results 
thus beg the question if properties such as the RKKY 
coupling, which are intrinsically linked to the magnetic 
properties of graphene, can accurately be described in 
a non-interacting electron picture. In this Letter, we 
will therefore investigate the effect of electron-electron 
interactions on the RKKY coupling, both in the bulk 
and in zigzag graphene nanoribbons (ZGNRs) where the 
zero-energy edge states 13| can significantly modify the 
RKKY behavior 0] . We will below show that even for 
small to moderate strengths of the electron interactions, 
the RKKY coupling in the bulk is qualitatively modi- 
fied and gets significantly more long-ranged than in the 
non-interacting electron picture. For ZGNRs the efi^ect 
is even more striking as any finite U causes the long- 
distance RKKY coupling to become distance indepen- 
dent. We thus conclude that it is imperative to include 



electron interactions when studying the RKKY coupling, 
and, by extension, any other properties closely related to 
the magnetic properties of graphene. 

More specifically, we will use the mean-field approxi- 
mation of the one-band Hubbard-JJ model for graphene 
and include magnetic impurity spins S = itSz which 
couples to a graphene atom with a Kondo coupling term 
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Here c^o- {c]^) annihilates (creates) an electron at site 
i with spin a, <i,j> means nearest neighbors, and 
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with aafj being the Pauli matrices, is 
the electron spin. The constants entering, apart from Jk 
which depends on the particular impurity moment, are 
the nearest neighbor hopping in graphene t = 2.5 eV and 
the on-site repulsion U. The value of U is hard to deter- 
mine exactly but, depending on the choice of exchange- 
correlation potential, U/t = 1 — 2 has been shown to be 
consistent with density functional theory (DFT) results 
[l^ . Below we are able to extract U/t — 2.1 when com- 
paring the RKKY coupling in a spin chain with DFT 
results [13]. The expectation value of the spin-resolved 
electron density rii^ = cJ^Cio- needs to be calculated self- 
consistently in Eq. ([T]) and gives the spin polarization 
density as sf = (n^-i- — ni^)/2. The Hubbard model 
has been employed before in the study of graphene and 
ZGNRs and has been shown to yield results consistent 
with first-principles DFT results [lilflTt. 

In standard RKKY perturbation theory [l^ the lead- 
ing interaction between two impurity moments at sites i 
and j is given by 
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with the effective RKKY coupling constant Jij propor- 
tional to the static spin susceptibility of the imbed- 
ding bulk. Here we will instead self-consistently solve 
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Eq. ([T]) for two impurity spins in a ferromagnetic (FM) 
and an antiferromagnetic (AFM) configuration, respec- 
tively and explicitly calculate the RKKY coupling as 
the energy difference between these two configurations: 
J^j = [E{FM) - E{AFM)]/2. More details of our 
method applied to non-interacting graphene can be found 
in Ref. [?] . 

Bulk impurities. — Figure [T] shows the magnitude of the 
RKKY coupling as function of impurity distance R along 
both the zigzag (a) and armchair directions (b) of the 
graphene lattice for several values of U /t. The RKKY 
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FIG. 1: (Color online) Dimensionless RKKY coupling 
\Jij\'^/Jk ^ function of impurity distance R in units of the 
lattice constant along zigzag (a) and armchair (b) directions 
for U/t = 0, 1, 1.5, 2, and 2.15 (increasing \ Jij\). A-A sublat- 
tice impurities (black, x) has FM coupling (Jij < 0) and A-B 
sublattice impurities (red, o) has AFM coupling {Jij > 0). 
Lines are only guides to the eye. 

coupling in the large _R-limit for non-interacting graphene 
is Jij oc [1 -t- cos(2kD ■ R)]/|R|3 with Jij < for A-A 
sublattice coupling, i.e. for impurities on the same sub- 
lattice, (black) and Jij > and three times larger for 
A-B (or different) sublattice coupling (red) Here 
ki) is the reciprocal vector for the Dirac points. Apart 
from minor effects due to a small R, these results are dis- 
played in the lowest black and red curves in Fig. [TJ The 
non-oscillatory i?-dependence for the armchair direction 
is a consequence of only sampling the cos-function at the 
graphene lattice sites. When including electron interac- 
tions these results are, however, qualitatively modified 
even for small U. For U/t — 1.5 (middle curve) essen- 
tially all evidence of the (1-1- cos)-oscillations is gone as 
is the factor of three difference between A-A and A-B 
sublattice coupling. Also, the power-law decay expo- 
nent a changes from 3 for U — to around 2.3 (2.6) for 
U/t = 1.5 and 1.9 (2.1) for U/t = 2 for the zigzag (arm- 
chair) direction. In fact, for U/t > 2 (uppermost curve), 
the armchair and zigzag RKKY couplings are equal and 
thus all lattice specific details have been washed out for 
such values of the electron interactions. With the mean- 
field quantum critical coupling for the AFM insulating 
state being Uc/t = 2.23 [19[ it is perhaps not surprising 



that the RKKY coupling becomes independent of the 
small length scale details close to this point. However, 
what is rather unexpected is that this "washing" out of 
the lattice details is clearly present even at such low val- 
ues as U/t = 1, a value which is very likely lower than the 
physical value of U in graphene. We thus conclude that 
including electron interactions is imperative when study- 
ing the RKKY interaction in graphene. Without them 
not only are the magnitude of the RKKY coupling grossly 
underestimated but, more importantly, the results do not 
even have a qualitatively correct i?-dependence. 

ZGNR impurities. — Within the non-interacting elec- 
tron picture we recently showed that for impurities along 
a zigzag graphene edge (A-A impurities) the RKKY in- 
teraction decays exponentially for large i?, but that, quite 
counterintuitively, smaller Jk gives a longer decay length 
Q. These results are a consequence of the extreme eas- 
iness by which an edge impurity can polarize the zero 
energy edge state. In contrast, for A-A impurities inside 
a narrow ZGNR, bulk properties of the RKKY coupling 
are largely regained, notably Jij oc J\/R? . The effect of 
the edge is thus only limited to edge impurities in the 
non-interacting limit. These results are shown in the two 
lowest curves in Fig. [2] for impurities along the edge (a) 
and inside the ribbon (c) for Jk = t (black) and Jk = t/10 
(red). When including electron interaction this picture 
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FIG. 2: (Color online) —Jij for Jk ~ t (black) and Jk = t/10 
(red) as function of impurity distance R for U/t — 0, 1, 1.5, 
and 2 (increasing — Jij) (a,c) and as function of U/t (b,d) for 
A-A edge impurities (a,b) and A-A impurities inside a narrow 
ZGNR of width W = 8/^/3a (c,d). Dashed lines shows Edw 
whereas dash dotted lines are equal to JkSsl with sf being 
the graphene polarization at the impurity sites but with the 
impurities absent. 

is dramatically changed. As is well established, any fi- 
nite U is going to spontaneously polarize the edge state 
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IgI 20j , and by extension the whole ribbon , thus mak- 
ing it harder for an impurity spin to influence the polar- 
ization of the graphene. The three upper curves in (a,c) 
are for U/t = 1, 1.5 and 2, respectively. As seen, the R- 
dependence completely disappears for R larger than a few 
unit cells for any physically relevant value of U and for 
all impurity sites in a narrow ZGNR. The i?- independent 
value of the RKKY coupling is analyzed as a function of 
U/t in Figs. [2i;b,d). As in the bulk, the FM impurity 
configuration is energetically favored for A-A impurities 
in ZGNRs, whereas the AFM configuration will require 
modification of the spontaneous graphene polarization to 
accommodate the impurity spin of the opposite orienta- 
tion to said polarization. There are two i?-independent 
limiting solutions for the AFM configuration of which 
the one with lowest energy will give an upper bound for 
the constant RKKY coupling. The first limiting solu- 
tion has a magnetic domain wall formed between the two 
AFM oriented impurity spins. The magnetic domain wall 
formation energy per edge Euw is equal to the RKKY 
coupling for this solution and its value, calculated within 
Eq. (jlj, is displayed with a dashed line in Figs. [5] (b,d). 
This limiting solution is not only independent of R but 
also of Jfc making it especially favorable at high -values 
which is also seen in Fig. [JJb). For smaller Jk it is, how- 
ever, more likely that the most favorable AFM solution 
is one where the impurity spins do not noticably change 
the polarization of the underlying graphene, not even di- 
rectly at the impurity site. The limiting AFM solution 
in this case is the unperturbed graphene plus the two 
impurities and has an energy 2JkSsf above that of the 
FM solution. Here sf is the graphene polarization at 
the site of the wrongly oriented impurity but in the ab- 
sence of impurities. This unperturbed limiting solution 
is also naturally _R-independent and its RKKY coupling 
is displayed with dash dotted lines in Figs. [2lb,d). Note 
that since sf is significantly lower inside a narrow ZGNR 
than on the edge, this solution yields a smaller RKKY 
coupling for impurities away from the edge. Also note 
that s| depends rather strongly on U, which causes both 
limiting solutions to increase sharply with an increasing 
U. The true RKKY coupling follow the lower of these two 
limiting solutions remarkably well for all four cases stud- 
ied in Fig. [21 including jumping form one to the other 
around U = 0.3t for edge impurities when Jk — t/10. 
The small deviations from the unperturbed limiting so- 
lution are due to some impact of the impurity spins on 
the graphene polarization which locally produces small 
changes in in favor of a lower total energy. The only 
real notable discrepancy is for edge impurities when Jk 
is large and U moderately small. Here the domain-wall 
limiting solution is not followed too closely but the sys- 
tem lowers its energy slightly by instead creating a local, 
half-circle shaped, domain wall around one of the wrongly 
oriented impurity spin. This solution naturally creates a 
spin- imbalance in the system as its domain wall does not 



propagate to other edge. Note that both of the limiting 
solutions described here are always present, and thus the 
qualitative RKKY behavior is the same, in any system 
which has a spontaneous polarization in the absence of 
impurities, such as is the case for ZGNRs for any finite 
U. For edge impurities we do not expect the width of the 
ribbon to change the RKKY behavior as both the spon- 
taneous edge polarization and Ejjw are weak functions 
of the ribbon width. However, for impurities inside a 
very wide ZGNR the spontaneous polarization inside the 
ribbon is going to be vanishingly small and bulk proper- 
ties should eventually be restored for wide ribbons and 
small Jk- We thus conclude that any finite electron in- 
teraction renders the long-distance RKKY coupling in a 
ZGNR _R- independent and linearly dependent on Jk for 
small Jk , but independent on Jk in the limit of large Jk ■ 
In addition, electron interactions make impurities inside 
a ZGNR behave similarly to edge impurities which is op- 
posite to the situation for non-interacting electrons. 

Determining U . — There exist some DFT results for the 
RKKY coupling in graphene 15|, |2]J but such studies 



are always limited to very small R unless chains (or lat- 
tices) of impurities are studied. Fig. [3] shows the RKKY 
coupling for A-A sublattice impurity chains along the 
zigzag direction separated a distance of 25 A as function 
of the impurity distance R along the chains. We see that 
for [/ = (lowest curve) characteristic non-interacting 
(l+cos)-type oscillations are present but the chain config- 
uration makes the RKKY coupling somewhat more long 
ranged than i?^"^. The oscillations however quickly dis- 




R (a) 



FIG. 3: (Color online) —Jij (black, x) for A-A sublattice 
impurities chains as function of impurity distance R along 
the chains for U/t = 0, 1, 1.5, 1.75, 2.1, and 2.15 (increasing 
magnitude). DFT results from Ref. [151 (dashed red, o) are 
scaled with a factor 0.5. Inset shows the power law decay 
exponent a as function of U /t (black, x ) with the exponent 
1.43 from Ref. [l5i] indicated with a dashed red line. 

appear and the decay exponent a decreases (see inset) 
with increasing U. DFT results using the hybrid func- 
tional B3LYP on the same chain structure is available in 
Ref. 15] and these results, scaled with an overall, unim- 
portant prefactor, are displayed with a dashed red line in 
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Fig. |3l There are no oscillations in the DFT results and 
the exponent a agrees when U jt = 2.1 which also yields 
very well matched results as indicative in the main plot. 
At a first glance, this might seem as a large value for 
the Coulomb repulsion in graphene since the mean-field 
AFM instability is at Uc/t — 2.23. However, one should 
keep in mind that multiple recent theoretical work have 
classified graphene as being very close to, if not even 
being, an insulator in vacuum due to strong Coulomb 
interactions l3"12|. Our results point to the fact that 
this state might be an AFM insulator which would be 
consistent with earlier results 12|. Our extracted value 



of U also agrees quantitatively with earlier estimations 
based on the B3LYP functional [31 ■ DFT calculations 
instead using the local density (LDA) or general gradient 
(GGA) approximations have yielded a somewhat smaller 
U/t - 0.9 - 1.3 [il It is well known that LDA 

suffers from electron self-interaction and therefore often 
underestimates U . B3LYP on the other hand explicitly 
contains an element of Fock exchange and thus tends 
to handle this deficiency better. This becomes especially 
important in strongly correlated systems but B3LYP can 
still reproduce the LDA results for weakly correlated ma- 
terials. 

With such high value of U it is also natural to ask 
about other possible electronically driven ordered states. 
?7c(AFM) increases with doping [l^ and thus undoped 
graphene is the strongest candidate for an AFM state. 
However, with increasing doping electronically driven d- 
wave superconductivity caused by spin-singlet nearest 
neighbor correlations appears for any Coulomb interac- 
tion [23|. Such correlations were already proposed by 
Pauling and others for the p7r-bonded planar organic 
molecules of which graphene is the infinite extension. 
With the Coulomb interaction extracted from the results 
in Fig. |3]one would need a chemical doping of = 1 eV 
to reach rc(SC) 5 K, a value which might be achieved 
with, for example, chemical doping. 

In summary we have shown that it is of vital im- 
portance to include electron-electron interactions when 
studying the RKKY coupling in graphene. Even rel- 
atively weak electron interactions qualitatively change 
the RKKY coupling to be significantly longer ranged 
and monotonically decaying in the bulk. In a ZGNR 
the change is even more pronounced and the R depen- 
dence entirely disappears. By comparing our mean-field 
results we have also been able to extract a surprisingly 
high value for the Coulomb interactions, demonstrating 
that graphene might be very close to an AFM insulat- 
ing instability. With such closeness to an AFM state it 
is rather natural that magnetic properties, such as the 
RKKY coupling, are going to be heavily influenced by 
electron interactions. 

The author thanks Sebastian Doniach, Jonas Fransson, 
Biplab Sanyal, Lars Nordstrom, and Eddy Ardonne for 
valuable discussions. 
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